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Focus_ Error = 0 
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Initialization 



304 

/ XJATindov/Center 
I YJ/vlndow Center 

X_Wldtri 

Y_Height 

Number_of_2_Slices 
t Z_Step_Size 
\ Start_Z value 



Move Z axis to Z_Position(N) = Start_Z value 
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Move Z axis to 
Z_PosJtion(N)- Z_Positlon(N) + 
Z_Step_Size 
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Acquire and 
Store Image(N) 
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Extract Image(N) 
focus Value F(N) 
Subroutine 
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Best_Focus = F(N) 
Best Z = N 
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Best_Z = 
Number_of_Z_Slices/2 
Focus_Error = 1 
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Set Bestjmag = Image(BesLZ) 
Set BestJrriage_Z = Z_Position(Best_Z) 
Seta Polnter_to_BestJmage 
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Extract Image(N) 
focus value F(N) 
Initialization 



PolnterJo_Stored Jmage ( 
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Convolve Image(N) with HorzSobel Filter to find 
Horizontal Edges creating lmage(N)_H 
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Convolve Image(N) with Vert_Sobel Filter to find 
Vertical Edges creating lmage(N)_V 
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SUM lmage(N)_H + lmage{N)_V to create 
lmage(N)_Focus 
Truncate values to range of 0 to 255 
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Calculate Variance F(N) in 
lmage(N)_Focus over the 
desired window 
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Initialization 
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Pointer_to_BestJmage 
Type_ ^Classification 
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Calculate Average_Gray_Background normalized for well! type by 
multiplying by Micro-well normalizer value 



Calculate Average_Gray_Classify_Window 
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Calculate Diffjmagejl 
If ABS[Pixel(x,y) - Pixel(x+Diff Sep,y+Diff_Sep)] > RagJThresh 
THEN 

Diff lmage_1_Pixel(x,y) = Set1 
ELSE 

DlffJmage_1_Pixel(x,y) = 0 
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Calculate Diff_lmage_2 
i If ABS[Pixel(x+DifLSep.y) - Pixel(x ( y+Diff_Sep)] > 
\ FlagJThresh 
THEN 

Diff Image 2_Pixel(x,y) = Set1 
ELSE 

DiffJmage_2_Pixel(x,y) = 0 
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Calculate Classifyjmage 
If Diff lmageJLPixel(x,y) = Set1 
OR 

DiffJrnage_2_Plxel(x,y) = Set1 
THEN 

Classify lmage_Pixel(x f y) ■ Set2 
ELSE 

ClassffyJmage_Pixel(x,y) = 0 
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Calculate the number_Pixels_Set2 
Calculate Total_Pixelsjn_Window 
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L_Gray = Average_Gray_Classify_Window/Average - _Gray_Background 
Score_Flag = Number_Pixels_Set2/TotaLPixelsjn_Window 
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CALL 

\Classify (Score Gray, Score _Flag) subroutine 
FIG 34a step 468 
Type_of_Classification as parameter 



X_BackGround_Center 
Y_BackGround Center 
X_BackGround_Width 
Y_BackGround_Height 
Micro-well type Identifier 
Micro-well normer value 
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RETURN 
to FIG 36 
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Classify 
(Score_Gray, 
Score_Flag) t 
.subrotine. 



470 



Initialization 
CLASS = 5 
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Calculate CLEAR 
If Score Flag< Clear_Flag_LT 
AND 

Score_Gray > Clear_Gray_GT 
THEN 
CLASS = 0 



Return to FIG 33 
Step 466 
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Calculate Light Precipitation 
If Score_Flag < Lgt_Precip_Flag_LT 
AND 

Score_Flag > Lgt_Precip_Flag_GT 
AND 

Score_Gray > Lgt Precip_Gray_GT 
THEN 
CLASS = 1 



Calculate Heavy Precipitation 
If Score_Flag < Heavy_Precip_Flag_LT 
AND 

Score_Flag > Heavy_Precip_Flag_GT 
AND 

Score Gray < Heavy Precip_Gray_LT 
THEN 
CLASS = 2 



Calculate Ugly Precipitation 
If Score Flag < Ugly_Precip J^agJJT 
AND 

Score_Flag > Ugly.Precip_Flag.GT 
AND 

Score_Gray < Ugly_Precip_Gray_LT 
THEN 
CLASS = 3 
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Calculate Micro_Crystals 
If Score Flag > Micro_Cry_Flag_GT 
AND 

Score Gray < Micro_Cry_Gray_LT 
. THEN 
CLASS = 8 



Calculate Crystal 
If Score_Flag < Crystal_Flag_LT 
AND 

Score_Flag > Crystal Flag_GT 
AND 

Score_Gray > Crystal_Gray_GT 
THEN 
CLASS = 9 
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Calculate GRAINY Precipation 
If CLASS = 8 
AND 

Score_Flag > Grainy_Flag_GT 
AND 

Score Gray > Grainy Gray_GT 
THEN 
CLASS = 7 
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\ Calculate a first set of Additional Image 

Features in Classifyjmage 
CALL Subroutine MVT_ToolsJ3lob_Analysis 
with pointer to Window in Classifyjmage 
Get 
num_found 
and for each found GET 
area(m), height(m), width(m), perimeter(m) 
location_X(m), location_Y(m) 
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m = 1 
num_spherolite = 0 
spherolite(m) = 0 
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Calculate Blob Ratios 
Circle_Like_HW(m) = height(m)/width(m) 
Circle_Like_AHW(m) = area(m)/(height(m)*width(m)) 
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Calculate Spherolite 
If CLASS = 8 OR9 
AND 

Circle_Like_HW(m) < Circle_Like HWJJpper 
AND 

Circle_Uke_HW(m) >Circle_Like_HW_Lower 
AND 

Circle_Like_AHW(m) < Circle Like_AHW_Upper 
AND 

Circle Like AHW(m) >Circle_Like AHW_Lower 
THEN 

num_spherolite = num_spherolite +1 
m = m+1 
spherolite(m) = 1 
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Calculate a second set Additional Image 
Features in Bestjmage 
Using the results from the previous 
MVT_Tools_Blob_Analysis on the 
Classify J mage 
m = 0 
num_phase_sep = 0 
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' ( Average_Gray_Background 
Pointer_To_BestJmage 
num_found 
area(m), height(m), 
width(m). perimeter(m) 
location_X(m) 
location_Y(m) 
num_spheroiite 
spherolite(m) 
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Usihg location_X(m), location.Y(m), 1/2(height(m)), 1/2(width(m) 
calculate the average gray value within this reduce window for each 
spherolite previously found and set phase_sep_Gray(m) = to the 

found value 
then Normalize it by calculating 
phase_sep_Gray(m) =phase_sep_Gtay(m) / 
.Average_Gray Background : 
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' phase_sep_Gray > phase_sep_Gray_GT 
THEN 
CLASS = 4 
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phase_sep_Gray_GT 
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Analyze previously measured features from FIG 34c to 
further separat CLASS 9 into sub-classes 9.0 through 
9.9 
set m = 0 

First test if CLASS = 9, if not then goto step 778 
if true then 
CLASS = 9.0 
if numjound =1 then goto step 768 
if not goto769 
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location_X(m) 
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HW_Ratio = Height(m) / Width 
(m) 

If HW_Ratio > needle_HW_GT 
or 

HW_Ratio < needle_HW_LT 
then CLASS = 9.2 
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If HW_Ratio > plate_HW_GT 
or 

HW_Ratio < plate_HW_LT 
then CLASS = 9.4 
THEN goto 773 ELSE goto step 778 




plate_HW_GT 
plate_HW_LT 
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Using location_X(m), location_Y(m), 1/2(height(m)), 1/2(width(m) ) 
calculate the average gray value within this reduce window for the Blob 
previously found and set Chunk_Gray(m) = to the found value 

then Normalize it by calculating 
Chunk_Gray(m) =Chunk_Gray(m) /Average_Gray_Background 
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\ If Chunk_Gray < Chunk_Gray_LT 
THEN 
CLASS = 9.6 
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If area(m) > Chunk_50_GT 
THEN 
CLASS = 9.8 
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If area(m) > gorgeous_GT 
THEN 
CLASS = 9.9 
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Return to FIG 33 
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For m = 0 to numjound 
Calculate for each blob 




.Average_Gray_Backg round 

Pointer_To_BestJmage 

numjound 

area(m), height(m), width(m), 

perimeter(m) 

location_X(m) 

location_Y(m) 

num_spherolite 

spherolite(m) 
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m = m + 1 
HW.Ratio = Height(m) / Width 
(m) 

If HW_Ratio > needle_HW_GT 
or 

HW_Ratio < needle_HW_LT 
then CLASS = 9.1 
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If HW_Ratio > plate J4WJ3T 
or 

HW_Ratio < plate_HW_LT 
then CLASS = 9.3 
THEN goto 798 ELSE goto step 812 
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plate_HW_GT 
plate.HW^LT 



NO 
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V ... 



Using location_X(m), location_Y(m), 1/2(height(m)), 1/2(width(m) ) 
calculate the average gray value within this reduce window for the Blob 
previously found and set Chunk_Gray(m) = to the found value 

then Normalize it by calculating 
Chunk_Gray(m) =Chunk_Gray(m) /Average_Gray_Backg round 



NO 
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If Chunk_Gray < Chunk Gray_LT 
THEN 
CLASS = 9.5 
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Chunk 



:_Gray_LT f 
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If area(m) > Chunk_50_GT 
THEN 
CLASS = 9.7 
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- | Chunk_50_GT 




FIG. 35b 



Return to FIG 33 
step 466 
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Type_of_Classrficati n / 
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Move to micro-well within 
selected micro-well plate 




-NO — i 
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acquire image and locate 




droplet 
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Move to droplet within 
micro-well in micro-well 
plate 



acquire droplet high- 
resolution image 
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CALL Specimen Auto Score & 
Classify Subroutine 
FIG 33 step 440 
Type_of_Classification as 
parameter 
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